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o 

^ , Abstract 

C^ . We perform a phase-space analysis of the cosmological 3-fluid problem consisting of a barotropic 

fluid with an equation-of-state parameter 7— 1, a pressureless dark matter fluid, plus a scalar field 
(j) (representing dark energy) coupled to exponential potential V — Voe\p{—KX(j)). Besides the 
^^ ' potential-kinetic-scaling solutions, which are not the unique late-time attractors whenever they ex- 

ist for A" > 37, we derive new attractors where both dark energy and dark matter coexist and the final 
O ' density is shared in a way independent of the value of 7 > 1 . The case of a pressureless barotropic fluid 

^Y' (7 = 1) has a one-parameter family of attractors where all components coexist. New one-parameter 

5h ' families of matter-dark matter saddle points and kinetic-matter repellers exist. We investigate the 

^^ stability of the ten critical points by linearization and/or Lyapunov's Theorems and a variant of the 

theorems formulated in this paper A solution with two transient periods of acceleration and two 
CN . transient periods of deceleration is derived. 

> 

o 

Q ! Observations have ever confirmed the so-called transient period of acceleration (TPA) |T|, in this regard 

cn ' any realistic cosmological model should include component(s) with negative pressure. In this work we 

consider a spatially flat Friedmann-Roberston- Walker (FRW) cosmology where the content of the uni- 
verse has two components with negative pressure and a pressureless component. This is the so-called 
S^ I 3-fluid problem consisting of a barotropic fluid with equation of state Py = {j— 1)Py where < 7 < 2, 

a pressureless dark matter density poM and a dark-energy-scalar-field <p coupled to exponential potential 
V{(j)) = Vbexp(— K'A^) with equation of state p^ = (O^Pij,, p^ = 0^/2 — ^(0) and p^ = 0^/2 + y(0). In 
this cosmological 3-fluid model — 1 < G)0 < 1; however, from a physical point of view, there is no com- 
pelling reason to constrain the values of ft)^ to the interval [—1, 1], in this regard it has been established 
that the teleparallel dark energy [2| cosmological model allows for ft)^ < — 1 |3|. 

Existing analytical methods [4] have failed [5] to produce exact solutions to the 3-fluid problem. 
Apart from some trivial solutions (power law inflationary solutions where the scale factor of the universe 
evolves as aoc t"" for m > 0), no exact solution to the 3-fluid problem seems to exist to our knowledge. 

In Ref. [6] we restricted ourselves to ordinary matter, where 1 < 7 < 2, and to positive exponential 
potentials and resorted to numerical approach by which we derived new solutions to the 3-fluid problem. 
The solutions where classified hyperbolic and trigonometric according to the value of A, this extends the 
classification made for the solutions to the cosmological 2-fluid problem |T| (consisting of a barotropic 
fluid plus a dark-energy-scalar-field <p coupled to exponential potential), which were first derived in \^. 
For the whole range of A , we were able to construct solutions with one TPA where the universe undergoes 
the deceleration-TPA-etemal deceleration expansions. No solutions with two or more TPA's were found. 
However, as we shall see later, solutions with one TPA and a late-time eternal acceleration expansion do 
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exist for the range < 7 < 2/3 (which were not reported in Ref. ||6]), where the universe undergoes the 
deceleration-TPA-TPD -eternal acceleration expansions (TPD: For transient period of deceleration). We 
shall also derive solutions with two TPA's and two TPD's for 7 approaching from below 2/3. Solutions 
with many TPA's and TPD's may exist too. 

Phase-plane and -space analyses of autonomous differential equations lead to specific solutions that 
may provide the late-time attractors or the early-time repellers. Both type of solutions are interesting and 
provide rich insight into the evolution of the universe. Prior to the determination of the exact solutions 
to the cosmological 2-fluid problem ||7l- IIT2I . phase-plane analyses of the 2-fluid problem were per- 
formed IIT3I - lITSl (for a general procedure see llT6l ) and led to the discovery of potential -kinetic-scaling 
solutions, which are the unique late-time attractors whenever they exist for A^ > 37. 

We shall carry a phase-space analysis of the autonomous differential equations governing the dynam- 
ics of the 3-fluid problem. Among the conclusions we reach (1) the stability of the scalar field dominated 
solution for A^ < min(3,37), (2) the stabihty of the potential-kinetic-matter scahng solution for 7 < 1, 
which are quantitatively different from the corresponding results for the 2-fiuid problem fT3l, and the ex- 
istence of (3) new attractors (the potential-kinetic-dark matter scaling solution and the potential-kinetic- 
matter-dark matter scaling solution) and (4) new repellers and saddle points. In Sect. [2] we derive the 
autonomous differential equations of the 3-fluid problem and some other useful formulas. In Sect. |3]we 
discuss and extend the methods used for the stability analysis. In Sect. |4] we derive the critical points, 
investigate their stability and and their cosmologic implications. We conclude in Sect. [5] 

2 Autonomous differential equations of the 3-fluid problem 

With similar notation as in fT3l, the dynamics of the three fluids in a spatially flat FRW universe, with 
a scale factor a{t) and a Hubble parameter H{t) = a/a, are governed by the very Eqs. (1) to (3) of [,13 J 
upon slightly modifying the first equation by adding the contribution attributable to dark matter 

^ = -y(Pr + Pr + PDM + 0^) = -y(7Pr + PDM + 0^), (1) 

-mypy, (2) 

(3) 

where F = AF /At. These equations are constrained by 

//2 = (fcV3)[p^ + PDM + (0V2)+V]. (4) 

From now on we consider only positive potentials V{<^) = Vbexp (— K"A0) where A > 0. The dimension- 
less variables 

^~ ~ y = ^=--.z = ^--, vv= "^^ , (5) 



Pr = 


-3H{Py + Py) 


= 


■ dV 



where F = HF' and F' = dF/dlna, reduce the system ([l}- ^ to the following system of three linearly 
independent autonomous differential equations 

x' = ^-Xy^ - 3x + \x[\ +x2 -/ + (^_ i),2], (6) 

y' = -^-Xxy + \y[l+x^ -y^ + {Y-\)z\ (7) 

z=-\yz + \z[\+x^-y^ + {y-l)z\ (8) 



where the variable w is solved by 



x^+y^ + z^ + w^ = l, (9) 



which follows from dUl. It is worth mentioning that the expression in the square parentheses in the 
system ^- ^ is positive or zero: 1 +x^ —y^ + (7— l)z^ = 2x^ + jz^ + w'^. To arrive at ^- dS) we used 

H'/H = -3{2x^ + yz^)/2. (10) 

The equation governing the motion of w is independent of A 

w' = 3w[x2-3;2 + (7-l)z2]/2. (11) 

In general (for all 7), the constraint ^ restricted the motion to within the unit solid 2-sphere of center 
at the origin: x^ +y^ + '^ <\. A necessary formula for the stability analysis is readily derived upon 
combining ^, ^ and (HJ setting a:^ +y^ + z^ = r^ 

(r2)' = 3(r2-l)(2x2 + 7z2_r2). (12) 

Another useful formula for the stability analysis and qualitative behavior of the solutions is derived 
upon eliminating the expression in the square parentheses in ([H) and (fTTI) 

z' w' 3,, , 
z w 2 

leading tco 

z2=l2wV(i-^), (13) 

where L > is a constant of integration. For a pressureless barotropic fluid (7 = 1), Eq. (fT2l) reduces to 
z^ = L^w^ which leads, using Q and setting I = L/\/Z?+T < 1, to 

x^+y^ + z^/f = l. (14) 

Thus for 7 = 1, the motion happens on an ellipsoid of revolution around the z axis in the phase space. 
The ellipsoid, which is inside the 2-sphere x^ +y^ + z^ < 1, does not contain all the trajectories for 7 = 1; 
as we shall see, there are some equilibrium points inside and outside the ellipsoid; there are also other 
trajectories corresponding to L = (=^ z = 0) and to L = 00 (=^ w = 0). 

The relative densities are defined by Q.^ = x^ +y^, Q.y = z^, O-tie = w^ and obey the conservation 
equation Q.^ +Q.j + Q.x)b = 1. Other relevant quantities are the parameter (O^ = {x^ —y^)/{x^ +y^) which 
is constrained by — 1 < ft)^ < 1 and the deceleration parameter q = —d/{aH^) which is, by the field 
equations, the same as —1 —H/H^ = — 1 —H'/H leading to 

lq = \+^[x^-y^ + {y-l)z^]. (15) 

Combining ( fT5l ) and (fT2l) we arrive at 

Q!j,^={2q-\)a^E (16) 

a; =[(2^-1) -3(7-1)]% (17) 

a; = (2<7-l)(l-a^) + 3(7-l)%. (18) 



'Eq. | |13I > is also derived upon combining Eqs. (18) and (19) of ID. 



3 The critical points (CP's) - Lyapunov's Stability and Instability Theo- 
rems (LST and LIT) 

The CP's are the equihbrium points (xc, yc, Zc, Wc) in the phase space obtained upon solving the nonlinear 
algebraic equations x' = 0,y' = 0, z' = 0, and w' = 0. To determine the stability of the CP's we proceed to 
the linearization of the system (O- dUl setting x = X +Xc,y = Y ->ryc, z = Z + Zc {w = W -irWc) where the 
new variables {X, Y, Z) still obey the full nonlinear system ^- ^. Upon linearization, the system Q- (H) 
is brought to the matrix form: 

{X',Y',zY=J,-{X,Y,Zf, (19) 

where {X, Y, ZY is the column matrix transpose of (X, 7, Z) and Jc is the 3x3 Jacobi matrix lITTll - |fT9]| 
evaluated at the CP {x^yc, Zc)'- 



Jc 



Jcu (-3xc + \/6A)jc 3(7-l)xcZc 

{3xc - \J\X)yc Jell 3(7- l)ycZc 

'iXcZc -3ycZc Jc33 



(20) 



where 27,ii = 3[- 1 + 3x2 -^2 + (7- y^^2^^ 27,22 = 3 - VGXx,. + 3x2 _ 9^2 + ^^_ y^^i^ 27,33 = 3[x2 - 

yl + ir-mzl-i)]. 

The test for stability of almost linear systems llTTI states that |[T7l[T9l if (1) all the eigenvalues of the 
nonsingular matrix 7, (det7c 7^ 0) have negative real parts, then the CP is asymptotically stable but if (2) 
any eigenvalue has positive real part, then the CP is unstable. If some eigenvalues are zero (det7, = 0) 
or have zero real parts (and still det7c / 0), we will employ appropriate arguments (among which Lya- 
punov's Theorems, LST ifTTl - |[T9l . and LIT ITFtI ) for the determination of the stability of the correspond- 
ing CP as the above-mentioned test is no longer valid ifTTlllQII . Mathematically speaking, we shall not 
make use of the notion of saddle points for a saddle CP is generically unstable. However, physically, we 
shall distinguish between a repeller and a saddle point. 

LST assumes the existence of a continuously differentiable function U{X,Y, Z) that is positive def- 
inite in a neighborhood & of the CP and has an isolated minimum at the CP, which is the origin in the 
new coordinates (X, F, Z): (Xcp, Yqv, Zqp) = (0, 0, 0). If further the derivative of U along a solution 
curve, U' = diUiX')' with / = 1 -> 3 and X^ = X, ^2 = y, X^ = Z, is negative definite on ^ (except at the 
origin): U'{X, 7, Z) < 0, then the CP is asymptotically stable. We are not concerned with the case where 
the CP is stable El- m. 

LIT ifTTll for 2-dimensional systems generalizes to higher dimensional systems in a straightforward 
way. It consists in finding a function U{X,Y, Z) that is continues on a domain & containing the CP, 
which is assumed to be isolated. The Theorem assumes that f/(CP) = and that there is at least a point 
Pq = (Xq, Yq, Zq) in each disc in &, of center CP, where U{Po) > 0. If U' is positive definite on & (except 
at the origin): U'{X, Y, Z) > 0, then the CP is unstable. 

The intuition behind LIT is as follows. Assume the above conditions are satisfied. In any disc Dg 
in & select a solution curve that starts al| Pq: X{N = 0) = Xq, F(A^ = 0) = Fq, Z{N = 0) = Zq. If the 
solution curve evolves from Pq to, say. Pi = (Xi, Fi, Zi) we must have L'^(7'i) > L'^(7'o) > since U is 
increasing along the solution curve. Now, since U{Po) > U{CP) = 0, this solution curve won't reach the 
CP in a finite or an infinite time N (otherwise U would decrease). Thus the CP is not asymptotically 
stable (an asymptotically stable point is a CP where any solution curve starting in its vicinity ends up at it 
as N ^ 00). Furthermore, the solution curve must leave the disc Dg since, otherwise, as N ^ 00, U ^ 00 
too, which is not possible as the continuity of U on Dg implies that it is bounded there. Thus, the CP is 



In a general problem, use t instead of A'. 



not stable; it must be unstable. LST works, in a sense, the other way around in that its hypotheses ensure 
that the solution curve approaches the CP as A^ — )• oo. 

In LIT, U need not be zero at the CP since one can add any positive or negative constant to U without 
modifying the condition of stability, and U{Pq) need not be positivqj: It is sufficient to have [/(CP) < 
U{Pq). a variant of LIT may be formulated as follows. If (1) U{x, y, z) is continues on a domain & 
containing the CP, which is assumed to be isolated, (2) in every disc centered at the CP [here the CP is 
not necessarily the origin of the coordinates (x, y, z)], there exists some point Pq = (xq, jo, Zo) such that 
[/(CP) > U{Pq), (3) U' is negative definite on ^, then the CP is unstable. 

In cosmology both LST and LIT are very useful. One is interested to stable CP's or attractors where 
different solution curves end up at regardless of their initial conditions. One is also interested to unstable 
solutions or repellers which represent starting or intermediate events. 

The main difficulty in applying LST and LIT is that there is no method how to find U . There are, 
however, some directions for that purpose ||T91 . However, the construction of U may be greatly simplified 
relying on the assumptions of LST and LIT concerning &. This will be illustrated in the following 
discussion. 

4 The critical points (CP's) - Stability analysis - Cosmological implica- 
tions 

We have counted ten CP's labeled from A to J. In the following we will provide the values of the CP's in 
the form (x^, yc, Zc, Wc), determine their stability conditions and discuss their cosmological implications. 
Following the classification made for the analytic solutions to the 2-fiuid problem [7 1, the solutions 
with A^ < 6 are called hyperbolic and those with A^ > 6 are called trigonometric. Due to different 
conventions, the value of A used in |6], Anum, is related to value of A used in this work by Anum = \/3A. 

A = (0, 0, 0, 1). For 7 7^ 1, Jc has at least one positive eigenvalue: {—3/2, 3/2, 3(1 — 7)/2}. This 
CP is unstable. For y = I, J^ is singular. However, it is straightforward to show that in this case the CP 
is also unstable. This is achieved upon linearizating ^ and (Q in which case we obtain the eigenvalues 
=F3/2 of opposite signs. 

Cosmologically, the only solution curves that may reach A emanate from B± with 3^ = and z = 0; 
otherwise, some solution curves (only those emanating from B+) may just get close to it, but do not cross 
it, as shown in Fig.[T] At this CP, all densities vanish for a dominant dark matter component Hde = 1, ft)^ 
is an indeterminate there and the universe undergoes a decelerating expansion with q = 1/2. 

B+ = (+1, 0, 0, 0), B= = (-1, 0, 0, 0). The matrix J^ has the eigenvalues {3, (6 - £\/6A)/2, 3(2 - 
7)/2} where £ = 1, —1 for B+, B , respectively, so they are generically unstable. They are also unstable 
in the special case 7 = 2 where Jc is singulaio as can be concluded from the linearization of ^ and (|7]i. 

These are the repellers, as shown in Fig. [Tjand Fig. |2] with a dominant kinetic energy, a decelerating 
expansion q = 2, and (O^ = I. For a steep potential, A > V^, B+ is a saddle point. 



■'The LST and LIT were firstly formulated to deal with the stability of autonomous differential equations where the CP is the 
origin of the new coordinates {X, Y, Z). This is no longer the case in other coordinate systems as (x, y, z). 

For 5+, if y = 2 and A = V6, we conclude to the instability upon applying LIT with U = aX and a > 0. The instability of 
the case 7=2 and any A can also be achieved considering ( I12t which becomes (r)'=3(r— l)(r^ — 2y ) . A solution curve that 
starts near the CP has r < 1 (the CP is on the sphere r=l and all solution curves are inside the sphere). Since y =Y <^r^ ^ 1, 
we have (r^)' < and thus the solution curve moves in the direction of decreasing r and never returns back to the CP where 
r = 1, which is then unstable. This is a first application of a variant of LIT which was formulated in the previous section. 



D = (X/VS, i/l-(AV6), 0, 0). This CP exits for A^ < 6 (the case A^ = 6 leads to the previous 
case). The eigenvalues of J^ are: {(A^ — 3Y)/2, X^ — 3, (A^ — 6)/2}. In the case det/c 7^ 0, the CP is 
asymptotically stable for A^ < min(3,37) and unstable for min(3,37) < A^ < 6 [this includes the cases 
(A^ = 3 and 7 < 1) and (A^ = 37 and 7 > 1)]. If det/c = 0, it is asymptotically stable in the cases (A^ = 3 
and 7 > 1) and (A^ = 37 and 7 < 1) upon linearizing [(jT) and dUl] and [Q and ^], respectively. 

There remains the case A^ = 3 and 7=1 where we expect the CP [in this case D = (1/a/2, 1/a/2, 0, 0)] 
to be asymptotically stable. We apply LST and select U of the form: U = a{X +Y)^ + {b — a)Y-^ + cZ-^, 
which is positive definite ifO<a<b and < c. The CP is an isolated minimum of U with U{CP) = 0. 
The directional derivative along the solution curves, U' = diU{X'y with / = 1 — )• 3 and X^ = X , X-^ = 
Y, X^ = Z, is evaluated using the r.h.s's of ©, ([7]) and dD after converting to new coordinates {X, Y, Z): 

U' = F{X,Y,Z) with F = -3{b-a)Y^ + 0[{x')^], (21) 

which is negative definite in the vicinity of the CP. We choose & to be any neighborhood of the CP 
(including the CP), where U' < 0, and not including other points of the surface S : F{X,Y,Z) = 0. This 
way we make U negative definitqj in &. With these choices we satisfy the hypotheses of LST, so the CP 
with the case A^ = 3 and 7 = 1 is asymptotically stable. 

For A^ < min(3,37), this CP is an attractor with a dominant scalar field component Q.^ = I, co^ = 
— 1 + A^/3 < and a deceleration expansion 2^ = 1 + A^. 

E = (0, 0, 1, 0). From the set of the eigenvalues of J,-, {-3(2 - 7)/2, 37/2, 3(7- 1)}, the CP is 
generically unstable if det/c 7^ 0. 

Now, consider the case where Jc is singular (det/c = 0). In the special case 7 = 2, the CP is unstable 
for the linearization of ^ and dD leads to Y' = 3Y , Z' = 3Z. The same conclusion is achieved from 
{X'^ + Y^+Z^y = 3{Y^ + Z^) + • • • > 0. For 7= 1 the CP is also unstable by LIT or upon hnearizing © 
and ^ which results in the eigenvalues ±3/2 of opposite signs. The case 7 = is stable since we have 

{x^+Y^+z^y = -3(y2 +z2) + • • • < 0. 

Thus, £■ is a matter dominant attractor (Q.y = 1) if 7 = 0, a saddle point if < 7 < 2, and a repeller 
if 7 = 2. With the parameter C0(j, remains undetermined, the state of the universe at E undergoes a de- 
celerating expansion if 2/3 < 7 < 2 or an accelerating expansion if < 7 < 2/3. This is a novel point 
because one may have a TPA without necessary having (at the same time) a minimum kinetic energy and 
a maximum potential energy as in the case of the 2-fluid problem ||6l[8l. In fact, at E both kinetic and 
potential energies are zero. Thus if, for < 7 < 2/3, a solution curve approaches the saddle point E then 
deviates to a CP (this would be the CP G), the TPA there (at E) may last longer than the TPA occurring 
away from saddle points. This is because a saddle point behaves partly as an attractor and partly as a 
repeller. This is in fact the case in Fig. |3]that is a plot, for A = \/6.3 and 7 = 0.6666, of twice the decel- 
eration parameter, 2q, and the kinetic and potential relative densities, x^ (dashed line) and y^ (continuous 
line). The parameter q crosses the N axis at: A^i = 4.82333, N2 = 7.984, A^3 = 10.3342, A^4 = 13.0676 
and N5 = 13.9537. This solution has thus two TPA's and two TPD's: The first and second TPAs are 
observed in the intervals Ni <N <N2 and A'^3 <N <N4, respectively, and the first and second TPD's are 
observed in the intervals N2 < N < N3 and N4 < N < N5, respectively. The first TPA starts at N = Ni, 
which is the moment where a:^ ~ and j^ ~ [x(A^i) = -0.075, 3'(A^i) = 0.106, z{Ni) = 0.991], that is 
the corresponding point on the solution curve is near E. The graph of 2q continues to oscillate for N > N5 
below the Une q = 0, this is a sign that solutions with many TPA's and TPD's may exist if careful choice 
of the parameters is carried out. Fig. S]is a similar plot with different inputs A = \/63 and 7 = 0.6. It is 



^This evokes the pendulum problem 1 18 1. Even if & were to include other points on the surface 5 (in which case U' would 
be negative semidefinite), we would still conclude that the CP is asymptotically stable (and not just stable). 



a solution with one TPA and one TPD which depicts a case with a TPA starting at the moment where x^ 
is minimum and y^ is maximum. 




Y 
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Figure 1 : Case A = 3, y = 1 .5. For these values of the parameters, H is the unique attractor. Solutions starting at S+ and E 
(for this value of 7, ii is a saddle point) get very close to A, which is a saddle point. All solutions starting in the vicinity of B± 
and E spiral to H. Those curves, which start in the vicinity of S± with j = and z = or in the vicinity of E with x = and 
y = 0, end up at A. Since A is unstable, any perturbations in the values of the coordinates cause the solution curve to continue 
its journey to H. 



F = (cos 6,0,sm6,0),0<6 <JC,y = 2. This unstable CP generahzes B± in that Q.y = sin^ d may 
assume any value between and 1; it also generalizes E. 

As Jc is singular, it is not possible to draw any conclusion concerning stability by linearization of the 
system ^- ^ for this CP. With 7=2, Eq. ^^ becomes (r^)' = 3(r^ - l){r'^ - 2y^). A solution curve 
that starts near the CP has r < 1 (the CP is on the sphere r = 1 and all solution curves are inside the 
sphere). Since y^ = Y^ <^ r^ Ki 1, we have (r^)' < and thus the solution curve moves in the direction 
of decreasing r and never returns back to the CP where r = 1, which is then unstable. This gives another 
application of a variant of LIT (see footnote IDl. 

Since 6 is not constrained, f is a new one-parameter family of kinetic-matter repellers. With a stiff 
equation of state 7 = 2, the initial density is shared between the barotropic fluid and dark energy, w^ = 1 , 
and q = 2 (decelerating expansion). Solution curves starting from F, which is represented by a semicircle 
in Fig. [51 reach H. 



G = {^/3/2(y/X), y/3/2y/ (2 -y)y/X, y/X^-Sy/X, 0). The corresponding solution for the 2-fluid 
problem 1 13] is a potential-kinetic scaling solution the stability of which does not depend on the value of 
7. For the 3-fluid problem we rather have a potential-kinetic-matter scaling solution, its stability depends 
on 7 as we shall see later soon. 

The eigenvalues depend on both (A, 7): (3/4){4(7— 1), 7—2 — a, 7 — 2 + a} where we define a = 
V^(2 - 7)[2472 + A2(2 - 97)]/A. This CP exists for V > 37 only. If det/c / 0, it is asymptotically stable 




Figure 2: Case A = 1.8, y = 1. For these values of the parameters, the vertical line J, through the end-points H = 
(V^/A, ^/TJl/X, 0) and G = {^/JJl/X, ^/fJl/X, \/l-(3/A2)), is the unique family of attractors {H and G are locally 
stable for A = 1.8, y = 1). / is the line through A = (0, 0, 0) and E = (0, 0, 1). There is a curve starting in the vicinity of 6_ 
which converges to a point on the line J. There are three curves starting in the vicinity of i?+. The upper and lower curves 
approach the line / then converge to different points on the line J. The intermediate curve, which corresponds to j = converges 
to the line /, which is a set of saddle points; any perturbation in the value of y causes this curve to end up at any point on the line 
J. [In this caption the coordinates of the CP's have been given on the form (xc, yc, Zc)-] 



for 7 < 1 and X^ >?>y. Furthermore, this CP is (a) a stable node for < 7 < 2/9 for alO A^ > 37 (the case 
7 = leads to E discussed above), (b) a stable node for 2/9 < 7 < 1 and 37 < A^ < 247^/(97 — 2), and 
(c) a stable spiral if 2/9 < 7 < 1 and A^ > 2472/(97-2). The CP is unstable for 1 < 7 < 2 and A^ > 37. 

If det7c = 0, this CP is asymptotically stable for 7 = 1 and A^ > 37 since the hnearization of (O 
and provides two negative eigenvalues: — (3/4)(l ibj8) with j8 = \/24 — TA^/A (a stable node for 
3 < A^ < 24/7 and a stable spiral for A^ > 24/7). For the remaining cases where Jc is singular (7 = 2 
or A^ = 37), the CP is unstable for 7 = 2 and A^ > 6 since near it we establish: {X^ + Y'^ + 2^)' = 

6(VA2 — 6Z + \/6X)2/A2 H > 0. For A^ = 37, the eigenvalues which result from the linearization 

of (O and d?]) are proportional to 7— 2 and 7— 1, ensuring asymptotic stability for 7 < 1 and instability 
for 1 < 7 < 2. For the case 7= 1 and A^ = 3, where det/f = 0, we recover the CP D = (l/\/2, l/\/2, 0, 0) 
which has been shown to be asymptotically stable. 

The fact that 7 < 1 — to ensure asymptotic stability of the CP — results in ft)^ = 7— 1 < 0, while in the 
case of the 2-fluid problem (O^, still given by the same formula, may have both signs. Q.^ and Q.y depend 
on both (A, 7): Q.^ = 'iy/X^, Q.j= \ — Q.^. With 2q = \ + 3(7— 1), the state of the universe approaching 
this stable point may undergo a decelerating expansion if 2/3 < 7 < 1 or an accelerating expansion if 
< 7 < 2/3. G is a saddle point for 1 < 7 < 2. 



H = (•^3/2/A, a/3/2/A, 0, \/A^ — 3/A). Here again the eigenvalues depend on both parameters 
(A, 7): (3/4){6(l - 7), -1 - j8, -1 +j8}. This CP exists for A^ > 3 only and it is asymptotically stable 
for 7 > 1 and A^ > 3. The CP is (a) a stable node if 3 < A^ < 24/7 (with 7 > 1) or (b) a stable spiral if 



This subcase exists also for the the 2-fluid problem but was not derived in 131 113] . In this subcase (0 < y < 2/9), A~ need 
not be smaller than 24y^/(9y- 2). 




2 4 6 8 10 



2q 
■0.0002 /f3 14" 



■0.0014 



N 



N 




2 2 

1 



0. 6 
0.2 




Figure 3: Case 1 



-12 



/63, r = 0.6666. For our initial conditions, at A' = 0, we took xq = -0.9999997, yo = 4.3569 x 10 

1— Xq— yjj. (Upper and lower left plots) Twice the deceleration parameter 2151. (Lower right plot) The kinetic and 

potential relative densities, x^ (dashed line) and y^ (continuous line). The parameter 2q crosses the N axis at: A'; = 4.82333, 
N2 = 7.984, N3 = 10.3342, N4 = 13.0676 and N5 = 13.9537. This solution has thus two TPA's and two TPD's: The first and 
second TPA's are observed in the intervals Ni < N < N2 and N^ < N < N4, respectively, and the first and second TPD's are 
observed in the intervals N2 < N < N^ and A'4 < N < N$, respectively. The first TPA starts at N = Ni, which is the moment 
where x ~0 and y ~ 0, that is the corresponding point on the solution curve is near E. 



X^ > 24/7 (with 7 > 1). If det/c = 0, it is asymptotically stable for 7 = 1 and A^ > 3 as the linearization 
of (l6ll and ^ leads to the same eigenvalues (3/4){— 1 — j8, —1 +j3}. For the case 7=1 and A^ = 3 we 
recover the CP D = (l/\/2, 1/V2, 0, 0) which has been shown to be asymptotically stable. 

The relevant parameters are Q.^ = 3/A^ which depends only on A, G)^ = 0, Udm = 1 — H^, and 
q=l/2. 

For the 2-fluid problem, G is the unique attractor for A^ > 37 (for all 7) llT3l . But G depends on 7, this 
means that the end-behavior of the solution depends on the nature of the barotropic fluid. We have seen 
that, for the 3-fluid problem, G is no longer stable for 1 < 7 < 2 but H is, which is a new attractor and 
does not depend on 7. Thus, no matter the barotropic fluid equation of state is, the universe's evolution 
ends up at the same state provided 1 < 7 < 2. As we shall se below, for 7 = 1 (A^ > 37), there is a line (a 
one-parameter family) of attractors all represented by the CP J. 



I = (0, 0, COS0, sind), <0 < Jt/2, y = 1. Jc is singular, however, the eigenvalues which result 
from the linearization of ^ and ^ are ±3/2, ensuring instability. 

This unstable CP generalizes E and A in that D-y may assume any value constrained by Hy + Q.um = 1 . 
with q = 1/2 and (O^ remains undetermined. Since 6 is not constrained, / is a new one-paiameter family 
of saddle points where only matter and dark matter are the nonvanishing components. / is represented by 
a vertical line in the phase diagram, which is the line AE of Fig. |2l 



/ = {y/3/2/X, ^yjJl/X, Zc, y/l-{3/X^)-zj), y = 1. The CP exists only for A^ > 3. With det7, = 
0, the CP is however asymptotically stable since the linearization of dSjl and ^ provides the two negative 
eigenvalues: — (3/4)(l it j8) (a stable node for 3 < A^ < 24/7 and a stable spiral for A^ > 24/7). 

Since is a free parameter, 7 is a new one-parameter family of attractors in which all components 





Figure 4: Case A = ^63, y = 0.6. For our initial conditions, at A'' = 0, we took xq = -0.999975, yo = 4.3569 x 10"'', 
zo = 0.000479. (Left plot) Twice the deceleration parameter 2q. (Right plot) The kinetic and potential relative densities, x 
(dashed line) andy (continuous line). The parameter 2q crosses the A^ axis at: A^i = 2.44678, A'2 = 2.94285 and A'3 = 6.14727. 
This solution has thus one TPA and one TPD: The TPA is observed in the interval Ni < N <N2 and the TPD is observed in the 
interval N2 < N < N^. The TPA starts at A'l, which is the moment where x^ is minimum andy^ is maximum. 

coexist: It is a potential-kinetic-matter-dark matter scaling solution where (0,j, =0, q = 1/2, ^0 = 3/A^, 
and n^ + ^DM = 1 — (3/A^). Q.^ and Q.om are both arbitrary and smaller than 1 — (3/A^). According 
to the last and most accurate observations [,20 J . D.^ = 0.721 ±0.015 (corresponding to A^ = 4.161), thus 
%< 0.279 ±0.015. 

In Fig. 121 / is any point on the line through A and E. Only solutions curves emanating from B± and / 
(including A and E) may reach J, which is represented by a vertical line in the phase diagram, this is the 
line //G of Fig. H 

5 Concluding remarks 

We have derived attractor solutions where both dark energy and dark matter coexist and the late-time 
density is shared according to Q.^ = 3/A^ and ^0 ±nDM = 1 in a way independent of the value of 7 > 1. 
The case of a pressureless barotropic fluid (7 = 1) is more interesting and has a one-parameter family of 
attractors where all components coexist with, as before, Q.^ = 3/A^ but ^.y + ^uM = 1 — (3/A^). New 
one-parameter families of matter-dark matter saddle points and kinetic-matter repellers were also derived. 
Solutions with one and two TPA's and one and two TPD's were derived and the quest for solutions with 
many TPA's is open. The ten CP's may be grouped into families as follows. 

(1) Repellers. These include B± and F and they are represented by the semicircle of Fig. |5] which 
includes £" if 7 = 2. Eqs. ^- ^ imply Q.'q^ = SQbe and Q.'y = 3(2 - Y)Q.y. Thus, for B± both relative 
densities, Ude and Q.y, increase at the beginning of the evolution against Q.(j, which starts decreasing. 
This applies to F too with the exception that D.y has a stationary value at the beginning of the evolution. 
In Fig. [J we plot three solution curves two of which come very close to A then converge to H. 

(2) Saddle points. If 7 = 1, these include all the points of the line through A and E (representing /). 
For < 7 < 2, £■ is a saddle point. They all, but E, have q = 1/2. 0)^ remains undetermined. For 7=1, 
as solution curves approach /, as shown in Fig.[2j all relative densities tend to become stationary as their 
derivatives vanish there by (IT6] l- (ITST l. Thus, / is a turning point. This is almost obvious from Fig.|2]where 
the two curves, which start from S+ and converge to two different values of J (here 7 is a one-parameter 
family of attractors which is a vertical line through H and G in the phase diagram), have their maximum 
values of z {^y ~ ^^) ^^ ^^^ vicinity of /. 

We have also noticed that, for < 7 < 2/3, a TPA occurs as the state of the universe approaches 
the intermediate state defined by the saddle point E (where both kinetic and potential energies are zero), 
which lasts longer than other TPA's occurring away from saddle points (where the kinetic energy has a 
minimum and the potential energy has a maximum). 
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Figure 5 : Case A- = 3, 7 = 2. For these values of the parameters, H is the unique attractor. The circle through 6+, E and B is 
the one-parameter family F of repellers plus B±. Any curve starting in the vicinity of this kinetic-matter repeller ends up at H. 

(3) Attractors. J and G form a set of attractors for 7 < 1 and X^ > Sy. / is a one-parameter family 
of new attractors represented by a vertical line in the phase diagram which extends from the point H = 
a/372/ A, a/372/A, 0) to the point G = {^/3j2/?i, ^/yi/X, y/l-{3/X^)). For 7 > 1, we have the 
potential-kinetic-dark matter scaling solution, H, which is a new attractor where the end-behavior of the 
universe's evolution does not depend on the barotropic fluid equation of state. G, the potential-kinetic- 
matter scaling solution, is stable for 7 < 1 but the universe approaching this late-time state undergoes a 
decelerating expansion, as it should be, only if (1 >)7 > 2/3. 

It is straightforward to see that the CP's correspond to power law solutions for the scale factor a{t), 
as is the case with the CP's of the 2-fluid problem HHIIll. If {Xc / and Zc / 0), we obtain from QO^ 
a{t) oc r with m = 2/ {6x1 + 'ijzD fT4l. 

For the case 7 = 1, it is interesting to give a qualitative description of the solution curves which 
lie on the ellipsoid (fT4l) . For 7=1, Eq. (fT2l) implies (r^)' = 3(r^ — l)(x^ —y^)- Thus, curves with 
higher kinetic energy densities (x^ > y^) move upward on the ellipsoid, in the direction of decreasing 
r i.e. decreasing Q.^ and increasing Q.y = L^Q.^yi, and those with lower kinetic energy densities (x^ < 
y^) move downward in the direction of increasing Q.^ and decreasing Hy = L^Hdm- The only critical 
point that lies on the ellipsoid is the point /ellipsoid = (3/2/A^, 3/2/A^, l\/\ — 3/A^), which also lies 
on the hne J through the points H and G, with Wc = a/1 —f-yjl — 3/A^ where < £ < 1 is still a 
free parameter, /ellipsoid lies on the segment of the ellipse that joins the points {\/^/2, l/\/2, 0) and 
(0, 0, 1). All solution curves end up, directly or spiraling, at /ellipsoid- Thus, for A^ > 24/7, since 7 is a 
stable spiral, the three relative densities, (^0, H^, Udm) undergo oscillations around their average values, 
(3/A2, £2(1 _ 2>/X'^), (1 - £2)(i - 3/A2)), respectively. 

This cosmological model of three fluids, consisting of a barotropic fluid with an equation-of-state 
parameter 7 — 1 , a pressureless dark matter fluid, plus a scalar field coupled to exponential potential 
V = Vbexp (— K"A0), offers more possibilities for alleviating the coincidence problem: The late-time state 
is a decelerating expansion if 7 > 2/3, (O^ < 0, and the late-time relative densities are constant (but depend 
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on A) or arbitrary with their values determined through observations only. 
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